Orbital Dynamics of Binary Boson Star Systems 
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We extend our previous studies of head-on collisions of boson stars by considering orbiting bi- 
nary boson stars. We concentrate on equal mass binaries and study the dynamical behavior of 
boson/boson and boson/antiboson pairs. We examine the gravitational wave output of these bi- 
naries and compare with other compact binaries. Such a comparison lets us probe the apparent 
simplicity observed in gravitational waves produced by black hole binary systems. In our system 
of interest however, there is an additional internal freedom which plays a significant role in the 
system's dynamics, namely the phase of each star. Our evolutions show rather simple behavior at 
early times, but large differences occur at late times for the various initial configurations. 



I. INTRODUCTION 

Boson stars are compact solutions to Einstein's grav- 
itational field equations coupled to a massive, complex 
scalar field. In addition to their interest as soliton-like 
solutions akin to so-called Q-balls 1], they also serve as 
a simplified model of an astrophysical compact object 
with which to probe non- linear regimes. 

In a recent work, we described an implementation of 
the Einstein equations that can evolve, among other pos- 
sible sources, boson star configurations 0. We stud- 
ied the head-on collision of boson stars with equal mass 
but possibly differing in phase and oscillation frequencies. 
The simulations revealed interesting effects arising from 
these possible differences. In particular, by adopting 
phase and frequency reflections the collisions displayed 
markedly different behavior, from merger to repulsion. 
In the present work we revisit these binaries in scenarios 
with non- vanishing angular momentum. 

Our interest in the present work is twofold. On the 
one hand, we want to understand the dynamics of the 
more complex scenario of orbiting stars. The presence of 
angular momentum in this case gives rise to richer phe- 
nomenology arising due to the orbiting behavior of the 
system and the peculiarities of single boson stars with 
angular momentum. As illustrated by stationary, ax- 
isymmetric, rotating boson star solutions whose angular 
momentum is quantized Q , rotating boson stars have an- 
gular momentum further constrained than their regular 
fluid counterparts. 

On the other hand, the system also probes the dy- 
namics of compact objects in general relativity. In this 
context it is particularly interesting to analyze to what 
extent rather simple analysis can shed light on the orbit- 
ing behavior and gravitational wave output of the sys- 
tem. We recall that binary black hole simulations have 
revealed, among other features, that the dynamics and 
consequent radiation from binary black holes can be rea- 
sonably well approximated by perturbative approaches 
well beyond the limited range of validity initially as- 
sumed 0, m. Such a result indicates that for binary 
black holes, the system undergoes a rather smooth transi- 



tion from a quasi-adiabatic inspiral regime to the merger 
and ring-down without significant nonlinearities becom- 
ing major factors. This rather simple behavior seems to 
be a consequence of the final object being another black 
hole. When this black hole forms, it encircles a larger 
region which 'cuts-away' portions of the spacetime with 
strong curvature and dynamics and the spacetime is de- 
scribed by a perturbed black hole. 

The binary boson stars considered in this work enable 
us to further examine compact binary systems in a differ- 
ent but certainly related context. For binary boson star 
cases, as for binary neutron star cases, the evolution at 
early stages demonstrates simple quasi-adiabatic behav- 
ior. However, as the stars approach each other, effects 
related to the stars' internal structure play a significant 
role. Here we study binary boson stars, analyze their 
dynamical behavior and extract the gravitational wave 
output for different cases. 

As we discuss, significant differences arise from two ef- 
fects. One effect is the same as that discussed in the head- 
on case arising from the effective interaction energy in- 
ducing significantly different gravitational potential wells 
when the stars come close to each other. 

The other effect concerns the possible end states of 
the merger given the presence of angular momentum. If 
the merger is to produce a single remnant star, presum- 
ably it would yield, asymptotically, a stationary, spinning 
boson star. However, for the system to settle into one 
of the spinning boson stars found in Q, it would have 
to shed angular momentum and/or energy to be consis- 
tent with one of the quantized states. There are several 
"channels" for the system to achieve this. One is through 
radiation of angular momentum and energy; this, how- 
ever, might require long dynamical times for the system 
to shed enough. In related orbiting systems, like binary 
black holes, only a fewpercent of the initial mass/angular 
momentum [1,0, S [1] is radiated. Thus, extraction of 
energy /angular momentum via gravitational waves from 
the merged object is likely a small effect here as well. 
Other channels, such as the break-up of the merged star 
and black hole formation, will play significant roles in the 
dynamics of the system as we illustrate in this work. 

For a particular configuration (i.e, masses and sepa- 
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ration of the stars), the evolution of boson/boson pair 
appears to be highly dependent on the value of the ini- 
tial angular momentum J^- This outcome contrasts with 
the results of our previous work ^2*1 which demonstrated 
that the (unboostcd) head-on collision (J^ = 0) results 
in a remnant black hole when the individual masses are 
large. 

The paper is organized as follows. In Section |TT] we 
briefly summarize the formalism for the Einstein-Klein- 
Gordon system, referring the reader to previous work for 
more details [2|. Section UTTl describes the procedure for 
setting the initial data. Section IIVI recalls some of the 
analysis quantities that are going to be used along this 
paper. Section |V] describes the numerical implementa- 
tion of the governing equations and the presentation of 
our results. We conclude in Section IVII with some final 
comments. 



obtained with the Galilean boost, similar behaviour is 
obtained when employing a Lorentzian one, since after 
some transient epoch the dynamics of the system is, for 
the most part, insensitive to the details of the initial 
data definition. The procedure for constructing the ini- 
tial data can be schematically represented in a few steps 

1. compute the initial data for the single boson star 
(i) centered at Xi, that can be described with the 
scalar field and the metric {</'^*\ ffi^ }• This step 
was explained in detail in 0. 

2. boost both the scalar field and the spacetime (and 
their derivatives) by performing the Galilean or 
Lorentz coordinate transformation. The boosted 
fields will be denoted by g^al^ ■ 

3. superposition of the independently boosted boson 
stars 



II. THE EINSTEIN-KLEIN-GORDON SYSTEM 

As described in the dynamics of a massive com- 
plex scalar field in a curved spacetime is described by 
the following Lagrangian density (in geometrical units 
G = c= 1) [l^ 



(1) 



Here R is the Ricci scalar, gab is the spacetime met- 
ric, 4> is the scalar field, (f) its complex conjugate, and 
^(I'/'P) = ™^ l<^P the interaction potential with m the 
mass of the bosonic particle which has inverse length 
units. Throughout this paper Roman letters at the be- 
ginning of the alphabet a, 6, c, .. denote spacetime indices 
ranging from to 3, while letters near the middle i, j,k, .. 
range from 1 to 3, denoting spatial indices. 

This Lagrangian gives rise to the equations determin- 
ing the evolution of the metric (Einstein equations) and 
those governing the scalar field behavior (Klein-Gordon 
Equations). We employ the same implementation de- 
scribed in to simulate our systems of interest, and the 
reader is referred to that paper for the complete details. 
The only difference is with the initial data described in 
the following section. 



III. INITIAL DATA 

We adopt a simple prescription to define the initial 
data describing two boson stars. Throughout this work 
this data will be obtained by a superposition of two stars 
with linear momentum. In order to define such data, we 
adopt the spacetime and the scalar field described in Q 
and boost each star either with a Galilean or Lorentz 
boost. While this data satisfy the constraints only ap- 
proximately, the constraint violation measured on the 
initial data is at or below the truncation error. Fur- 
thermore, although the runs presented here have been 



(j, = ,^(i)(f -xl) + 0(2)(-_3,-) 

gab 



(2) 



where g^'"*^ = diag(— 1, 1, 1, 1) is the flat metric in Carte- 
sian coordinates. The field (j)^'"'^ (the unboosted scalar 
field) is dictated by the type of boson star considered, 
that is, either a boson star or an anti-boson star. The 
prototype of the single boson star {i) used throughout 
this work has M^*) = 0.5 and radius r'^I = 24M(*) 
(defined as the radius containing 95% of the mass) , 
so its compactness is comparable to a soft neutron star 
(M(*V^^*^ - 0.04). We also note that these boson stars 
are field configurations of a single complex scalar field. 



IV. ANALYSIS QUANTITIES 

Throughout this work we define the center of each star 
as the location of a local maximum of the energy den- 
sity p = n° n'' Tab for each temporal slice. Due to the 
U(l) symmetry of the Lagrangian density (H]), there is a 
conserved Noether current defined by 



r = -2.9"" 9,0] . (4) 

The conserved Noether charge iV, associated with the 
number of bosonic particles, can be expressed as 



(5) 



The ADM mass and the angular momentum of the 
system are computed as 



M = ^— hm 

167r r^oo 



g'' [djg^k - dkg^j] J^" dS (6) 
gj^ tvK] dS{7) 
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where stands here for the unit outward normal to the 
sphere. 

The gravitational radiation is described asymptotically 
by the Newman-Penrose ^'4 scalar. To analyze the struc- 
ture of the radiated waveforms it is convenient to decom- 
pose the signal into —2 spin weighted spherical harmonics 
as 

M r *4 = ^ C/,„r'li,^ (8) 

where the factor r is included to better capture the l/r 
leading order behavior of ^'4. 

We also focus on integral quantities that are indepen- 
dent of the specific basis of the spherical harmonics, such 
as the radiated energy and the radiated angular momen- 
tum. Using both the decomposition ([8]) and the orthonor- 
malization of the spherical harmonics, these expressions 
can be written as 

l.rri 

^ = -3^E™(MA,™Wi?;,„.w]), (10) 

where we have used the adimensional quantities 
1 /■* 

Di,m = mJ Ci^t') dt' (11) 
1 /■* 

El,rn = J A,m(i') dt' . (12) 

V. SIMULATIONS AND RESULTS 

Our starting point is an equal-mass binary system with 
initial velocities and separation chosen so that a sim- 
ple Newtonian analysis would predict the system to be 
bound. While both stars have equal mass, we exploit 
the freedom in setting the sign of the frequency and the 
phase of the stars to concentrate primarily on two differ- 
ent cases. The first case consists of two identical boson 
stars (the BB case), and the other consists of a boson 
star paired with its antiboson partner (the BaB case). 
One could also consider a boson star interacting with a 
copy of itself offset in phase, but we defer this to future 
work. Here we begin by surveying the dynamical behav- 
ior of the BB case as its angular momentum is varied. 
This allows us to examine the possible phenomenology 
of the system. Next we concentrate on two particular 
cases with total initial orbital angular momentum (ori- 
ented along the z axis) given either by Jz = 0.8M^ or 
Jz = l.lIvP. These values have been chosen so that the 
angular momentum is either below or just above the low- 
est allowed quantized value of the (rotating) stationary 
axisymmetric boson star, which obeys 

Jz = nN (13) 



with n an integer. Roughly speaking, if A^^*) denotes the 
Noether charge for a stable star (where i ~ 1 denotes 
the first star and i — 2 the other), then its mass is ap- 
proximately given by Af^'^ ~ A^'^^-'/m, being = ±1 
whether it is a boson or an antiboson jtar and m the 
boson mass (see for instance Fig. 3 of [3]). As a result, 
the Noether charge evaluated on the initial hypersurface 
is approximately 

N= — ^-(eiM^'^+e2M(^A. (14) 
m m \ J 

Here we consider m = 1 and M'^^^ = M^^^ = 0.5, so 
that the initial Noether charge for the boson/boson case 
is ~ Af^ while for the boson/antiboson pair we have 
precisely A^ = 0, which can be exploited to analyze the 
possible outcomes of a merger (or close interaction) of 
the stars. 

In the BaB case, since the total Noether charge must 
be zero, the interaction can not give rise to a single sta- 
tionary regular object, so the final stage can not be the 
spinning boson star described in Q . A more intuitive ex- 
planation relies on observing that a prospective remnant 
boson star would no more likely have positive charge than 
negative and so cannot exist. Therefore, either annihila- 
tion, dispersion or the (unlikely) production of multiple 
pairs of boson/antiboson stars is to be expected in the 
cases where a black hole is not formed. 

In the BB case, the Noether charge is non-zero, which 
by itself does not rule out the possibility of a single reg- 
ular object as the final outcome. However, angular mo- 
mentum considerations do shed further light on this is- 
sue. Notice that the boson/boson binary cases we con- 
sider have orbital angular momentum with values slightly 
above and below the Noether charge. In the latter case, 
since the total angular momentum is below the minimum 
allowed value, an outcome describing a single star with a 
mass comparable to the total mass would seem impossi- 
ble. In the former case, while forming a single boson star 
is possible, the system must shed enough angular mo- 
mentum/energy to satisfy the relation l|13p . As we will 
see below, the dynamics of the system reveals a a rather 
involved behavior. We arrange our presentation so as to 
discuss the case with smaller angular momentum first, 
and then concentrate on the higher momentum case. 

We study the dynamics of these two cases using fi- 
nite differences to approximate the equations of motion 
within a distributed and adaptive infrastructure. For 
the simulations of the angular momentum survey in Sec- 
tion IV Al the computational domain covers the region 

g [— 200Af, 200A/] while the stars begin within about 
40Af of each other, where Af ~ 1 is the total ADM mass. 
This separation, which is larger than what we adopt in 
the following subsections, allows for a larger angular mo- 
mentum in a bounded binary system and thus our sur- 
vey is more exhaustive. Within this domain, finer grids 
are placed dynamically according to an estimate of trun- 
cation error provided by a self-shadow hierarchy. Typ- 
ically, refined regions contain both stars and the span 
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between them with even finer regions tracking the indi- 
vidual stars providing a minimum grid spacing for each 
star of Aa; = 0.125M and a maximum grid spacing far 
from the stars of Ax = 4.0M. 

For the other simulations, the computational domain 
extends to a;* G [-400M, 400M] although the stars begin 
within about 32M of each other. Such a large coarse grid 
allows us to extract gravitational radiation far from the 
dynamics. We compute surface integrals (e.g. the ADM 
mass) and ^'4 at extraction surfaces located at r^xt — 
150M, 180M and 210M, where the grid spacing is given 
always by Aa; = 2M. 
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FIG. 1: Phenomenology for the boson/boson interaction with 
angular momentum (specific to the family of initial data de- 
scribed in the text) . There are two main difi'erent behaviours 
in the bound case, the formation of either a black hole or 
a rotating bar. The rotating bar, depending on the angular 
momentum, will either split into two objects or it will dis- 
perse on a longer time scale. Animations obtained from the 
simulations for some of these cases are available online uM. 



A. Survey results 

In this subsection we describe briefly the dynamics 
found in the binary boson star system as we vary the ini- 
tial velocities of the stars. Here, we make a few general 
remarks about these results, but we should note that we 
have yet to fully explore the parameter space. However 
in the limited region so far explored, a rich phenomenol- 
ogy is uncovered by our simulations. We choose stars 
with masses AfW = 0.5 near the maximum mass allowed 
on the stable branch M„iax — 0.63. Further, while the 
stars begin centered at (a;, y, z) = (0, ±16M, 0), their ini- 
tial speeds are perpendicular to their relative position 
vector, {vx,Vy,v^) = (ivboost, 0, 0). 

We find that for low speeds, the stars approach and 
form a spinning black hole as one would expect for two 
bound massive compact objects which, in isolation, are 
already close to the unstable branch. For stars with very 
large initial speeds, the system is unbound and the stars 
continue to fly away from each other. 

The interest lies in the intermediate range of speeds 
for which the stars interact nonlinearly. In this regime, 
as one increases the initial velocities the stars merge into 
a short-lived rotating object which finally splits into two 
pieces that are ejected. Continuing to increase the initial 
speeds, one finds the ejecta decreasing in size with a con- 
comitant increase in dispersion of scalar field. Eventually 
as one increases the velocities one observes the formation 
of a central black hole. 

The various transitions in parameter space just de- 
scribed are sketched in Fig. [I] (and representative ani- 
mations can be found online [ill). Such phenomenology 
demands more thorough study. However, we defer to a 
future work this more detailed investigation as the asso- 
ciated computational costs are very high. In the present 
work our focus is to obtain the gravitational wave sig- 
nature for representative cases and compare the results 
with those of other compact orbiting systems. 



B. Small Angular Momentum (J^ = 0.8M^) 
Boson/Boson pair 

We consider a binary boson/boson star system with 
initial angular momentum of = 0.8 M^. As the evolu- 
tion proceeds, the stars approach, orbit about each other 
for about half an orbit before merging into a single object. 
This object however, after spinning as a single entity for 
about t ~ 200M, splits into two identical objects which 
fiy apart from each other. Thus, this case corresponds 
to what we have labeled as the spinning bar/scattering 
region of Fig. [TJ 

The trajectories are plotted in Fig. [21 where the coordi- 
nate position of the stars' centers are shown for different 
times. A careful inspection of the scattered stars indi- 
cate their mass is half what the initial stars have. Fur- 
thermore, a Newtonian calculation would indicate that 
the final speeds of the stars are above the escape speed if 
no other gravitating source were present. Most of the 
remaining mass was dispersed during the merger and 
spinning, and, although a small fraction crosses the ex- 
traction surfaces, most of the scalar field is found in an 
extended halo surrounding the central merger region. 

Therefore, the system can be regarded as undergoing a 
non-trivial scattering in which the initial stars give rise to 
an unstable merged object. The instability of this object 
likely results from having insufficient angular momentum 
to settle into even the lowest allowed rotating, stationary 
boson star. As a result the smaller stars are 'kicked-out' 
in a configuration whose angular momentum is slightly 
below the initial one. The details of this scattering can 
be seen in Fig. [3l where several snapshots of the energy 
density at different times are plotted. The remnant scalar 
field remains in an extended halo around the origin. 

The emitted waveforms show a dominant Z = 2, to = 2 
mode in the spin-weighted decomposition of Mr\E'4 com- 
puted on the surface extraction at r = 210M, as demon- 
strated in Fig. [31 Its qualitative features are reminiscent 
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FIG. 2: Boson/boson pair (J^ = O.SM^). The position of the 
center of the star, centered initially at {x,y) = (0, ±16M), 
during the merge and the scattering. The stars form a single 
rotating object for T — 200M before splitting in two equal 
objects that are ejected towards infinity. 




of the waveforms produced in zoom- whirl orbits of com- 
pact objects[T3|. The flux of energy and angular momen- 
tum are also plotted in Fig. 31 

Because the stars are rapidly moving towards the ex- 
traction surfaces, and errors are propagating from the 
outer boundary, the waveforms become untrustworthy 
after t = 800. Nevertheless the dominant portion of the 
gravitational wave output in the system is clearly seen. 
Additionally, as indicated in the figure, the amount of an- 
gular momentum radiated is significantly stronger than 
that of the radiated energy, although both of them rep- 
resent a small fraction of the initial quantities. 

Additionally, this dominant mode and the energy flux 
are computed at different extraction radii in Fig. [5] (top), 
showing the expected wave-like behaviour of the Weyl 
scalar and a convergent behaviour of the energy flux 
as the extraction radius is more distant. The effects of 
having the surface extraction too close to the stars are 
shown in the same figure (bottom), where the energy flux 
is plotted for the three different radius. The tail of this 
quantity, much more sensitive than to distance effects, 
converges to zero as the extraction surface is placed far- 
ther away. 

In our previous work, we presented convergence tests 
in which the resolution was increased by a fixed amount. 
These tests indicated the code converged as expected. 
Here, as a complementary test, we examine the solu- 
tion's behavior when we vary the refinement criterion. 
In particular, when the estimate of the truncation error 
exceeds a user defined threshold e, refinement is added 
there. Notice that this test is not designed to examine 
the convergence rate of the solution, but rather to study 
the implementation's ability to adjust the grid structure 
so that the solution's error stays below a given one. By 
decreasing the value of e the grid structure effectively in- 
creases the resolution where needed. Fig. [5] illustrates 
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FIG. 3: Boson/boson pair (Jz — 0.8M^ ). Snapshots at the 
z — plane of the energy density p for different times. As the 
stars come closer and merge, the maximum value of p grows 
significantly. The boson stars keep together for t — 200Af and 
then they fly away with signiflcantly less mass than initially. 
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FIG. 4: Boson/boson pair fJz = 0.8M^;. Different modes 
of the ^4 (normalized by M r) at the extraction radius r — 
210M as a function of time at the top and the flux of energy 
and angular momentum at the bottom. Notice that the flux 
of angular momentum is about two orders of magnitude larger 
than the flux of energy. 
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FIG. 5: Boson/boson pair (Jz = 0.8M^). The dominant 
mode C2,2 of Mr^4, and the flux of energy dE/dt extracted 
at spheres with different radii. The plots have been shifted 
horizontally by the radial distance between the spheres. 



the integral of the energy density and the principal mode 
of the 5*4 decomposition for three different values of e. 
As is evident from these graphs, the quality of the ob- 
tained solution improves as the threshold is decreased 
and displays convergent behavior. 



Boson/ Antiboson Pair 

In contrast to the BB case, the total Noether charge 
for the boson star and its antiboson partner is zero. One 
implication of this cancellation is that the system can- 
not settle into a single remnant star. In spite of this 
preclusion, we observe the stars merging into some sort 
of asymmetric, inhomogeneous rotating object. 

Starting from an initial configuration identical to the 
one from the previous subsection with the exception that 
one star is an antiboson star, the behavior is illustrated in 
Fig. [7l Snapshots of the energy density show the stars or- 
biting very close while much of the scalar field disperses 
from the origin. The system settles into a larger inho- 
mogeneous object, apparently a rotating composite of a 



boson and antiboson star. The trajectories of these com- 
ponents are plotted in Fig. [51 

The radius of this object, computed from the en- 
ergy density, varies with time. For instance, at about 
t ~ lOOOAf its size has grown to about five times the size 
it had shortly after the merger, Lq- The expansion con- 
tinues to about t ~ 1500M when its size is about QLq. At 
this point the expansion halts and turns around, decreas- 
ing its size to approximately SLq- Afterwards this size 
stays roughly the same. Most of the scalar field energy 
density does not disperse away but remains contained in 
a big sphere centered at the origin. Although we can 
see that the energy density crossing the surface extrac- 
tion is about two orders of magnitude larger than in the 
BB pair with the same angular momentum, the escaped 
scalar field is still a small portion of the total one. 

As before, we can compute 5*4 and decompose it into 
spin-weighted spherical harmonics. In Fig. [5] the lowest 
modes are plotted, extracted at Text — 150Af , showing 
again a clear Z = 2,m = 2 dominant mode. The energy 
and angular momentum fluxes are also plotted in the 
same figure, where it is clear that the pattern is different 
from the BB pair in Fig. 21 Additionally, once again, the 
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FIG. 6: Boson/boson pair (J^ = O.SAf^ j. The integral of the 
energy density Mp and the dominant mode C2.2 of Mr'^4, (ex- 
tracted at r = 210M and shifted horizontally hy t = 210M) 
for three different truncation error thresholds e. 



flux of angular momentum is two orders of magnitude 
larger than the energy flux. 



C. Large Angular Momentum (J^ — l.lAf^) 
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Here we consider both the case of two boson stars or- 
biting each other (the BB case) and a boson star and 
its antiboson partner (the BAB case) for large initial an- 
gular momentum. In contrast with the previous case, 
here the BB pair has angular momentum which exceeds 
the Noether charge, thus a possible end state is a sin- 
gle boson star with angular momentum. Alternatively, a 
spinning black hole could be produced. Since Jz — 
this is a delicate outcome because this would be close to 
an extremal Kerr black hole, unless a significant amount 
of angular momentum is shed. 

The situation with the BaB case, while less involved, 
still has interesting features. Since the Noether charge in 
this case is zero, a single boson star cannot be produced. 
Thus one expects that either a black hole forms or the 
scalar field disperses in such a way that the total angular 



FIG. 7: Boson/antiboson pair (Jz = 0.8M^). Snapshots at 
the z — plane of the energy density p for different times. 
Notice that the maximum of p decreases while the radius of 
the star increase. 



momentum, up to the typically small radiated amount, 
is accounted for. As we will see, the simulations reveal 
the resulting outcome is the formation of a black hole. 

To illustrate the dynamical motion of the stars, we plot 
the coordinate location of one of the stars in Fig. [TO] for 
the two cases. The boxes/circles indicate the location of 
the star in the BB/BaB cases at the same interval of an 
asymptotic observer's time. Notice that while early on 
the stars behave roughly in the same manner, as they get 
closer the BB case 'speeds' up and merges earlier than 
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FIG. 8: Boson/antiboson pair (,h = 0.8M^). The position 
of the extremes of the Noether density, centered initially at 
{x,y) — (0, ±16Af), which define roughly the location of the 
boson and the antiboson. 



the BaB. This can be understood by our analysis indi- 
cating the effective Newtonian potential for the BB case 
is deeper than that of the BaB case (see the Appendix 
ofi). 

The dominant mode Z = 2, to = 2 of the spin- weighted 
decomposition of M is shown in Fig. [11] for the BB 
pair and Fig. [T^] for the BaB pair. Notice that the final 
speeds are around v = 0.15 c, that is roughly twice the 
initial speed. Since the BB pair orbits faster before merg- 
ing, the associated gravitational wave output is stronger 
than the BaB case. The energy and angular momentum 
fluxes are shown only for the BaB pair because the fluxes 
for the BB pair were very similar. 

Following , we can compare three different measures 
of the orbital frequencies obtained from the numerical 
simulations, which can be defined under some assump- 
tions. The first one is the coordinate frequency Wc, ob- 
tained from the numerical trajectories of Fig. [TOl This 
quantity is susceptible to gauge effects because it depends 
on the dynamical coordinates used during the evolution. 
The second one, ujd, is extracted from the dominant 
mode of ^4, which can be estimated by 



ujd = Im 

TO 



(15) 



The third one is computed by assuming a Newtonian 
binary in a circular orbit. In this case, the standard 
quadrupole formula gives the approximate I = 2 mode of 
the inspiral waveform 



C2,±2(t) =32 



where /i = M'^^^ M'^'^^ / M is the reduced mass and 4>{t) 
the accumulated phase of the orbit. Then, the frequency 
^NQC computed from the Newtonian quadrupole circular 
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FIG. 9: Boson/antiboson pair (J^ = 0.8M^;. The principal 
modes of at the further extraction radius (r = 210M) 

as a function of time at the top and the flux of energy and an- 
gular momentum at the bottom. The reflections and bound- 
ary effects become dominant around t — lOOOM, prevent us 
to characterize the radiation of the final inhomogeneous spin- 
ning object. 



orbit approximation can be computed from (fTB|) as 



M UJNQcit) 



M 
32 n 



3/8 



\C2,±2{t)\ 



(17) 



These three measures of the orbital frequency are 
shown in Fig. 1131 where there is a pretty good agree- 
ment between the coordinate frequencies tOc and lu]j. The 
frequency lunqc shows also a remarkable agreement in 
the first stages, taking into account that the initial data 
is only near the Newtonian quasi-circular orbit configu- 
ration. However, at the last part of the plot, the fre- 
quency ujnqc differs significantly from the others be- 
cause its approximation does not consider interactions 
other than purely gravitational, namely scalar field in- 
teractions. The same quantities are shown for the bo- 
son/antiboson pair in Fig. 1141 
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FIG. 10: Boson/boson and Boson/antiboson pairs (J^ — 
I.IM^ ). The coordinate position of the Galilean-boosted bo- 
son stars for the BB and the BAB pairs. The stars complete 
roughly one orbit before merging. 
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FIG. 11: Boson/boson pair (^J^ = I.IM^). The real part of 
the dominant mode 6*2,2 of M r for the BB pair, com- 
puted from the numerical waveform at two different extrac- 
tion surfaces and from the numerical trajectories, by using 
the quadrupole formula. 



VI. COMPARISONS AND CONCLUSION 

In this work we have studied the different phenomenol- 
ogy arising from the addition of angular momentum in 
the boson star binaries, going further than in our pre- 
vious work 0. First, a survey of the possible different 
behaviours, -for equal-mass binaries- was performed by 
varying the initial angular speed of the stars. This sur- 
vey revealed the significantly different outcome that can 
arise due to the restrictions single stationary boson stars 
have as far as their allowed angular momentum. Next, 
several cases were studied in detail for the BB and BaB 
pairs, paying special attention to configurations resulting 
in orbits. Here the radiation of the system was computed 
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FIG. 12: Boson/antiboson pair (J^ — I.IM^). The real part 
of the dominant mode C2,2 of M r ^'^ for the BaB pair at 
the top, computed from the numerical waveform at two dif- 
ferent extraction surfaces and from the numerical trajecto- 
ries, by using the quadrupole formula. The flux of energy 
and angular momentum are at the bottom, extracted at the 
outer extraction surface and compared with the result from 
the quadrupole formula. Notice that the flux of angular mo- 
mentum is at least two orders of magnitude larger than the 
flux of energy. 



and few salient features compared with binary black hole 
simulations. A detailed comparison between binary black 
holes, binary neutron star and binary boson stars is un- 
dergoing and will be presented elsewhere [l^. 

It is worth stressing an important consequence of the 
survey performed. Except for some special cases, almost 
all initially "bounded" configurations (as dictated by a 
Newtonian analysis) give rise to an evolution describing 
either the formation of a black hole or the dispersal of 
the scalar field. Thus it would appear quite difficult to 
form a single, rotating boson star from the collision of 
high mass (and high compaction ratio) boson stars. Of 
course, one could try forming rotating boson stars with 
lower mass boson stars, or with other forms for the inter- 
action potential (e.g. terms quartic in the field modulus). 
However, our suspicion is that such formation, if it occurs 
at all, will turn out to be very much non-generic. 

The different particular cases studied showed rich fea- 
tures. The low angular momentum cases end up with 
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FIG. 13: Boson/boson pair (J^ = I.IPI^). The orbital fre- 
quency computed three different ways; ujc from the numerical 
simulation, lud from the dominant mode of ^'4, and ujnqc 
assuming the Newtonian quadrupole circular orbit approxi- 
mation. 
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FIG. 14: Boson/antiboson pair (J, = l.lA/2/ The same 
quantities as in Fig. ll3l but for the boson/antiboson pair. The 
lower values of the frequencies at the final times, with respect 
to the boson/boson pair, is due to the weaker effective inter- 
action between the boson/antiboson stars. 



a dispersed scalar field and a low-mass remnant object 
at the center. Even more interesting was the realization 
that in the BaB case, a closely rotating pair consisting of 
a boson/antiboson is induced. 

For the large angular momentum case, both configu- 
rations (BB and BaB) evolve to a spinning black hole 
after dispersing/radiating part of its mass and angular 
momentum. We computed the gravitational waveforms 
produced by the merger for both cases and compared 
them with the waveform computed using the quadrupole 
formula. Interestingly, while the stars are quite close ini- 
tially and only perform about an orbit before merging, 
the computed waveforms show significant agreement with 
those obtained from the quadrupole formula. Moreover, 
even in the more violent merger of the BB pair, there is 
a smooth transition between the spiral and the merger, 
as has already been observed in binary black hole merg- 
ers. Hence this system represents another example of 
compact binary coalescence for which the nonlinearities 
of the merger are barely manifest in the gravitational 
waveforms. 
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